Observation of B° D^r]' and S° D*^r]' 
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We report the observation of B*^ — > D^r]' and the first observation of B^ — » D'^ri' , using 140 
fb~^ of data collected at the T{4S) resonance with the Belle detector at the KEKB asymmetric 
energy e+e" collider. We find the branching fractions to be B(B° D'^v') = [1-14 ± 0.20 (stat) 
-o.vj (syst)]xlO"* and B(B° D*°i^') = [1.21 ± 0.34 (stat) ±0.22 (syst)] x 10"" with significances 
including systematic uncertainty of 8.9 and 5.3 standard deviations, respectively. 



The colour-suppressed decays £)(*)o /jO have 

been observed for /i" = T:^,ri,uj and 0, S 0> and 
recently also 77' 4J . The rates are all larger 

than originally expected in generalized factorization pj. 
The large branching fractions may be explained by fi- 
nal state rescattering ^ or non-factorizable diagrams in 
perturbative QCD [7|. Hard-coUinear effective theory 8] 
cannot explain the rate for D(*)°7r°; other decay 

modes have not been investigated in this framework. 

Both the colour-allowed decay B^ D^*^+h^ 9] and 
the colour-suppressed decay B^ D'-*^^h^ proceed via 
the emission of a virtual W~ . In the former case the 
W~ can decay without colour restrictions, however, in 
the latter case the d quark from the W~ decay has to 
carry colour that matches the colour of the d quark orig- 
inating from the . This allows colour-singlet formation 
only in one out of three cases resulting in suppression by 
a factor of i (therefore colour-suppressed) . The tree level 
diagrams for both decays are shown in Fig. ^ The con- 
tribution of the VF-exchange diagram is usually assumed 
to be neghgible j^]- 
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FIG. 1: Colour-allowed (a) and colour-suppressed (b) tree 
level Feynman diagrams for B D'^*^h. 

The branching fraction for B^ D'^rj' and an up- 
per limit for the B^ — > D*^r]' decay have recently been 
published by BaBar 0. In this Letter we present the 
first observation of B^ D*'^r]' and an observation of 
B^ D^rj' with more than 5S statistical significance. 



This analysis is based on a 140 fb ^ data sample con- 
taining 152 X 10^ BB events collected with the Belle de- 
tector at the KEKB ^ e+e" colhder. KEKB is an 
asymmetric energy collider (3.5 GeV on 8 GeV) that 
operates at the T(4S') resonance (y^ = 10.58 GeV) 
with peak luminosity of nearly 1.4 x 10^'^'* cm~^s~^. 
B^B^ and B^B^ pairs are assumed to be produced with 
equal rates. For signal and background simulation the 
QQ [13 Monte Carlo (MC) generator was used and the 
GEANT3 package was used for detector simulation. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer consisting of a three-layer silicon vertex detector 
(SVD), a 50-layer central drift chamber (CDC), an array 
of aerogel threshold Cerenkov counters (ACC), a barrel- 
like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of 
CsI(Tl) crystals (ECL) located inside a superconducting 
solenoid coil that provides a 1.5 T magnetic field. An iron 
flux-return located outside of the coil is instrumented to 
detect Kl mesons and to identify muons (KLM). The 
detector is described in detail elsewhere [l^. 

Charged tracks with impact parameters less than 0.5 
cm radially and less than 2 cm in z (the z axis is anti- 
parallel to the positron beam direction) with respect to 
the interaction point and with momenta larger than 100 
MeV/c are selected. Kaon and pion mass hypotheses are 
assigned to charged tracks on the basis of a likelihood 
Cx/tt that is obtained by combining information from the 
CDC {dE/dx), ACC and TOF systems. The likelihood 
ratio TZk = CK/{^Tr + Ck) ranges between (pion-like) 
and 1 (kaon-like) and we require TZk > 0.6 for kaon and 
TZk < 0.4 for pion candidates originating from the heavy 
meson yielding an efficiency of about 88% in both cases 
and TZk < 0.9 for pion candidates from the light meson 
side corresponding to an efficiency of 99%. In addition, 
we reject tracks that are consistent with an electron or 
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muon hypothesis. 

Photon candidates are selected with a minimum en- 
ergy requirement of 50 MeV. Neutral pion candidates 
are reconstructed by combining two photons with invari- 
ant mass between 115 MeV/c^ and 152 MeV/c^, which 
corresponds to ±2.5 standard deviations in terms of the 
mass resolution. In addition, we require the momentum 
of the tt" to be greater than 400 MeV/c in the laboratory 
frame. 

The 7] meson is reconstructed by combining two pho- 
tons. An asymmetric 77 mass window was chosen, corre- 
sponding to ^3 standard deviations. The mass of the 77 
meson is then constrained to its nominal value from the 
Particle Data Group (PDG) Q- 

We select p candidates from two oppositely charged 
pion candidate tracks in a mass region of —220 MeV/c^ 

< M(7r+7r-) - m(p) < 170 MeV/c^, where rup is the 
nominal mass of p° mesons from the PDG and M(7r+7r^) 
is the mass of the p candidate. In addition, we require the 
transverse momentum of the daughter pions, Pt{tt), to 
be greater than 300 MeV/c. This requirement suppresses 
around 40% of the continuum background while retaining 
86% of the signal candidates. 

The 77' meson is reconstructed in its two dominant de- 
cay channels 77' — > rjir'^Tr' and 77' pj. For both subde- 
cays the 77' candidates are required to be within a mass 
window of 0.94 GeV/c^ < M(77') < 0.975 GeV/c^ and 
have a center-of-mass (CM) momentum of Pc a/ ('7') > 1-7 
GeV/c. The pions from the rj' decay must have a trans- 
verse momentum in the labratory frame of Pt{t^) > 0.1 
GeV/c. For the latter subdecay, the photon from the 77' 
is required to have an energy greater than 0.2 GeV in the 
lab frame. This selection alone suppresses around 85% 
of background while retaining 85% of signal candidates. 
A weak requirement on the helicity of the 77 meson of 
1/1(77)1 < 0.97 is applied, where ^1(77) is the cosine of the 
angle between the 77' momentum and the direction of one 
of the decay photons in the 77 rest frame. An 77' mass 
constraint is applied for the final fitting results. 

The mesons are reconstructed in three different de- 
cay channels, D° K-Tr+, D° R-tt+tt^ and £)° 
if^7r+7r^7r+ with mode-dependent ±3(T mass windows. 
The — > i^^TT+Tr" decay channel has an additional re- 
quirement on the tt" momentum of f*(7r°) > 0.4 GeV/c 
in the lab frame to suppress continuum background. A 
mass constraint yields a significant improvement to 
the momentum resolution for the — > K^tt^tt^ decay 
and is applied for all decay channels for consistency. 

The D*^ mesons are reconstructed combining a 
meson with a tt" or a 7. We require the mass dif- 
ference of the D" and D*^ to satisfy 0.136 GeV/c^ 

< M{D°Tr°)-M{D") < 0.148 GeV/c^ and 0.131 GeV/c^ 

< M{D°-f) - M{D°) < 0.156 GeV/c^ In addition, we 
require the momentum of the D'^ in the lab frame to be 
greater than 0.8 GeV/c. 

The B'^ mesons are reconstructed combining an 77' me- 



son and a D° or meson. Two kinematic variables are 
used to extract the meson signal: the energy differ- 
ence AE = Eb — fi'bcam and the beam constrained mass 
^^bc = V^bcam " -^1 ' whcrc -Bbeam IS the beam energy 
and Eb and Pb are the reconstructed energy and momen- 
tum of the candidate in the T(4S') rest frame. The 
events that satisfy the requirements Mbc > 5.2 GeV/c^ 
and \A.E\ < 0.25 GeV are selected for further analysis. 

After these selection criteria, the two major back- 
ground sources are continuum e+e^ qq (where q = 
u, d, s, c) and b ^ c decays. 

Several event shape variables (defined in the CM 
frame) are used to distinguish the more spherical BB 
topology from the jet-like qq continuum events. The 
thrust angle 9t is defined as the angle between the pri- 
mary 5° decay daughters and the thrust axis formed by 
all tracks not from the same 5° meson. Jet-like events 
tend to peak near |cos0t| = 1, while spherical events 
have a fiat distribution. The requirement | cos 6*^1 < 0.9 
is applied prior to any other event topology selections. 

Additional continuum suppression is obtained by us- 
ing modified Fox- Wolfram moments |0| and the angle 
Ob between the fiight direction of the reconstructed B^ 
candidate and the beam axis. A Fisher discriminant 
(jF) is formed by a linear combination of cos 6t, 
S± and five modified Fox- Wolfram moments. S± is the 
ratio of the scalar sum of the transverse momenta of 
all tracks outside a 45° cone around the 77' direction 
to the scalar sum of their total momenta. Probabil- 
ity density functions (PDFs) are obtained from signal 
and background MC data samples. These variables are 
then combined to form a topological likelihood function 
Cc =PDFc(cos 6'b)-PDFc(.F) where c = signal (s) or con- 
tinuum background (qq). Signal follows a 1 — cos^ {9b) 
distribution while continuum background is uniformly 
distributed in cos9b- We select signal-like events by re- 
quiring a high likelihood ratio TZc ~ Cs/{Cs + Cqq) to 
suppress continuum background. For channels with an 
77' ^"7 decay an additional variable cos9t-c, which is 
the angle between the 77' momentum and the direction of 
one of the decay pions in the p rest frame, is included for 
better signal-background separation. The TZc require- 
ments are optimized using signal and continuum back- 
ground Monte Carlo samples and are found to be strongly 
mode dependent and ranging from 0.1 to 0.925, where in 
general more restrictive constraints are applied for decays 
with 77' — !■ p^j and looser constraints for decays including 
a.D*°. _ 

For events with multiple candidates, the best can- 
didate is selected on the basis of the for the vertex fit 
to the charged pions from the 77'. If multiple D*^*^ meson 
candidates remain after this selection, the best candidate 
is selected using — iM{D^*^) — 77i23(.) )^/(t^, where 
M{D^*y) is the candidate mass, m]j{,} is the nom- 
inal mass and a is the resolution of reconstructed 
mass. Multiple candidates in the heavy meson side 
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appear for less than 2% of the events for the — > D'^rj' 
decay but in the order of 20% for the i?" D*^r]' decay. 
Backgrounds from other B decay modes such as B' 



£)(*)t 



P~ 



and B — s- D^*^°a-^ are studied. These de- 



cays are suppressed by vetoing events where the p or 
ai could be constructed from the B° candidate daugh- 
ters with -0.1 GeV< AE < 0.08 GeV and Mbc > 5.27 
GeV/c^. The vetos suppress about 10% -30% of the BB 
background after the final selection criteria, while retain- 
ing over 98% of the signal MC events. For decays with 
r]' — > rj'K^TT~ , which have less background contamina- 
tion, the remaining BB contributions are modeled with 
a single two-dimensional smooth function, obtained from 
a large MC sample. For decays with -q' p^^ we divide 
these backgrounds into two groups, the B^ D^*^^a'^ 
decays and all others, both are again modeled with two- 
dimensional smoothed histograms. 

The signal PDFs are represented by a Gaussian plus a 
bifurcated Gaussian (a Gaussian with different width on 
either side of the mean) in AE and a bifurcated Gaus- 
sian in Afbc- AH PDFs are extracted from Monte Carlo 
simulations. For B'^ D^t]' events the feeddown from 
— > D*^!]' is modeled using 2-dimensional smoothed 
histograms. Continuum events are represented by a first 
or second order polynomial in AE and an empirical back- 
ground function introduced by ARGUS ^3 for Mhc- Sig- 
nals should peak around Ai? ~ GeV and Mbc = 5.28 
GeV/c^. Correction factors accounting for the difference 
between MC and data, determined from a study of the 
high statistics decay mode B" v'^s, are applied to 
the mean and width values of the signal shapes for both 
AE and Mbc. 

The reconstruction efficiencies are determined from 
a large sample of signal MC events, and range from 
0.6% — 11.2% for the 18 different decay channels, where 
the subdecay branching fractions are not included. Cor- 
rection factors due to differences between data and MC 
are applied for the charged track identification, photon, 
TT^ and 77 reconstruction. 

The signal yields (Ns) are extracted using extended 
unbinned maximum-likelihood fits simultaneously per- 
formed in AE and Mhc- An extended likelihood function 
is: 



LiNs,NB^) = 



NsPs{x^)+J2NB^PB^ix,) 



(1) 



where N is the total number of events, i is and in- 
dex running over the events, Xi IS cL vector of the AE 
and Mbc values for each event and Ps and Pb^ arc the 
probability density functions for signal and background, 
respectively, and the index j runs over all background 
sources. The signal yield Ns and background contents 
Nb are determined by maximizing the L(Ns, Nb ) func- 
tion in the {Ns, Nb ) manifold, where the Nb defines a 



j-dimensional submanifold of all different backgrounds. 
The statistical significance of the signal is calculated as 
E = y^2hi(Xfn^i^(7A))j where imax and Lq denote the 
maximum likelihood value and the likelihood value at 
zero branching fraction, respectively. 

The signal and background normalisations are floated 
in the fit while other PDF parameters are fixed to val- 
ues determined from MC studies. When combining two 
or more modes we fit for a common branching fraction 
instead of the signal yield for each mode. In the case of 
B^ D^T]', the contributions from the Z?*° feeddown are 
fixed to the branching fraction obtained in this analysis 
and the number of D*^ events is calculated with a mode 
dependent "feeddown-efficiency" , which is derived from 
MC simulations of B'^ D*'^ri' decays reconstructed as 
B^ — ^ D^rj' . The branching fractions together with the 
yields, the reconstruction efficiency for each decay and 
the significances are listed in Table HI The AE and Mbc 
projections of the fits are shown in Fig. 12 



TABLE I: Signal yields, efficiencies ^. ei5(-Df'), with 
the index i running over the D*-*' subdecay modes), branching 
fractions B and significances E. 



Decay mode 


Yield 


Efficiency 


B 


E 






(xlO-^) 


(xlO-^) 






49.4 ±8.7 


28.3 ±0.8 


1.14 ± 0.20 lo;i3 


8.9 




24.3 ± 6.8 


13.0 ± 0.4 


1.21 ±0.34 ±0.22 


5.3 



The major sources contributing to the systematic un- 
certainty of the branching fraction measurements are the 
systematics of the correction factors for the AE and Mbc 
shape parameters (1.9% - 3.4%, depending on the de- 
cay mode considered), the choice of the two-dimensional 
smooth functions (4.1% - 12.5%, this is estimated by fit- 
ting the data sample with one-dimensional functions in- 
stead, half the deviation is taken as the error, this system- 
atic error partly includes the systematics from the PDFs) 
and the uncertainty of the reconstruction efficiency of 
charged tracks (1.0% per track). Other systematic er- 
ror contributions are the uncertainty of the shape of the 
PDFs, which is quantified by varying each parameter of 
the PDFs by ila of its nominal value. The changes in 
the yield are added in quadrature, resulting in an error 
of 0.9% - 1.8%. The uncertainties in the reconstruction 
efficiencies for all 77, tt" and photons together give a sys- 
tematic error in the range of 2.9% - 3.8%. The particle 
identification fake rate and efficiency uncertainties for all 
charged tracks yield a systematic error of up to 0.4%. 
The efficiency is estimated from MC data samples with 
an accuracy of 2.5% - 5.0%. By requiring different TZc 
selection criteria the systematic error from this source is 
evaluated to be 2.0% - 8.9%. The subdecay branching 
fraction uncertainties as listed by the PDG result in 
an error of 1.5% - 3.3% and the number of BB events 
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FIG. 2: B° ^ D'^ri' decay /^E (a) and Afbc (b) projection 
plots and _B" — > D'^-q' decay /S.E (c) and M^c (d) projection 
plots. Points with error bars represent the data, the solid 
line is the signal contribution, the dash-dotted line is BB and 
the dotted line is continuum background. The two shaded 
areas represent the D*" D-k'^ (dashed line, filled area) and 
j~f*o _j jj^ (dotted line, hatched area) feeddown. The dashed 
line is the sum of all contributions. The projection plots are 
for signal regions of all other variables, the signal regions for 
and A/bc are: -0.055 GeV < < 0.05 GeV and 5.27 
GeV < < 5.29 GeV. 



that are recorded with the Belle detector is known to 
an accuracy of 1%. The systematic error contributions 
described above are added in quadrature and result in 
a total systematic error of ^if 5% for D^rj' and 

for W ^ D*^i^'. 

In summary, we observe the decays D'^i]' and 

— !■ D*^ri' with significances including systematic un- 
certainty of 8.9 and 5.3 standard deviations, respectively. 
The latter is observed for the first time. Our branching 
fraction for D^rj' is about one standard deviation 

below another recent measurement 0. Both branching 
fractions are higher than early theoretical predictions but 
are in agreement with recent work 0. 
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